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Abstract 

We show that the evolution of the number density of rich clusters of galaxies breaks the 
degeneracy between 0 (the mass density ratio of the universe) and cr 8 (the normalization of 
the power spectrum), a 8 fl 0 - 5 ~ 0.5, that follows from the observed present-day abundance 
of rich clusters. The evolution of high-mass (Coma-like) clusters is strong in H = 1, low-a 8 
models (such as the standard biased CDM model with a 8 ~ 0.5), where the number density of 
clusters decreases by a factor of ~ 10 3 from z = 0 to z ~ 0.5; the same clusters show only mild 
evolution in low-fl, high-cry models, where the decrease is a factor of ~ 10. This diagnostic 
provides a most powerful constraint on Q. Using observations of clusters to z ~ 0.5 — 1, 
we find only mild evolution in the observed cluster abundance. We find = 0.3 ± 0.1 and 
cr 8 = 0.85 ± 0.15 (for A = 0 models; for fl + A = 1 models, = 0.34 ± 0.13). These results 
imply, if confirmed by future surveys, that we live in a low-density, low-bias universe. 

subject headings : galaxies : clusters - galaxies : evolution - galaxies : formation - 
cosmology : theory - cosmology : observation - dark matter 


1. Introduction 

The observed present-day abundance of rich clusters of galaxies places one of the 
strongest constraints on cosmology (Bahcall and Cen 1992, White et al. 1993, Eke et al. 
1996, Viana and Liddle 1996, Pen 1996): a 8 H 0 5 ~ 0.5 ± 0.05, where <x 8 is the normalization 
of the power spectrum on 8 h _1 Mpc scale (reflecting the rms mass fluctuations on this 
scale), and is the present value of the cosmological density parameter. This constraint is 
degenerate in U — cr 8 ; models with 0 = 1 and <t 8 ~ 0.5 are indistinguishable from models 
with VI ~ 0.25 and cr 8 ~ 1. (A cr 8 ~ 1 universe implies no bias in the distribution of mass 
versus light, since cr 8 (gal) ~ 1 is observed for galaxies; a a 8 ~ 0.5 universe, on the other 
hand, is highly biased, with mass distributed more widely than light). 

In the present paper, we show that a study of the evolution of the number density of 
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rich, massive clusters as a function of redshift will break the degeneracy between O and a 8 
and determine each parameter independently. The growth of high mass clusters depends 
strongly on the cosmology — mainly O and a 8 (e.g., Press and Schechter 1974, Peebles 1993, 
Cen and Ostriker 1994a, Jing and Fang 1994, Eke et al. 1996, Viana and Licldle 1996). In 
low-density models, density fluctuations evolve and freeze out at early times, thus producing 
only little evolution at recent times (2 1). In an Id = 1 universe, the fluctuations start 

growing only recently thereby producing strong evolution in recent times: a large increase 
in the number density of massive clusters is expected from z ~ 1 to z = 0. The evolution 
is so strong that finding even ~ 1 — 2 Coma-like mass clusters at z ~ 0.5 over 10 3 deg 2 
of sky would rule out an 11 = 1, a 8 — 0.5 model, where only ~ 10 -2 such clusters would be 
expected (§3). 

We investigate in this paper the evolution of the mass function (MF) of clusters (the 
number density of clusters above a given mass) for various cosmologies using large-scale 
N-body simulations, and compare the results with cluster observations to 2 ~ 0.5 — 1. 


2. Model Simulations 

We investigate the evolution of the cluster MF in five cosmological models using large- 
scale numerical simulations. The models include: Standard Cold Dark Matter (SCDM; 
ff = 1), normalized to a present-day mass fluctuation on 8h -1 Mpc scale of a 8 = 1.05 
(consistent with the COBE microwave background fluctuations on large scales); a biased 
SCDM (ff = 1), with low normalization (cr 8 = 0.53), which fits the present day cluster 
abundance (but is inconsistent with the COBE normalization); a low-density, A dominated 
CDM model (LCDM), and an open CDM model (OCDM); and a mixed, hot and cold 
dark matter model (MDM, 0=1, 0^=0.3; from Cen and Ostriker 1994b, normalized to 
crg=0.6, consistent with the present day cluster adundance). Table 1 summarizes the model 
parameters. All models, except for the COBE-normalized SCDM (<jg — 1), are consistent 
with the present day cluster abundance (Bahcall and Cen 1992, White et al. 1993, Cen 
and Ostriker 1994b, Eke et al. 1996, Pen 1996). The SCDM a 8 ~ 1 model over-produces 
the number of massive clusters by an order of magnitude. All models except for the biased 
SCDM (<jg ~ 0.5) are also consistent with the COBE normalization (Bunn & White 1996). 

A large-scale particle-mesh code with box size 400 h -1 Mpc was used to simulate the 
evolution of the dark matter in the models. A large simulation box is needed in order to 
produce a significant number of the rich but rare clusters (<; 10 -5 h 3 clusters Mpc -3 ). The 
simulation box contains 720 3 cells and 240 3 = 10 71 dark matter particles, with a particle 
mass of 1.3 x 10 12 ffh -1 M Q . In each simulation, clusters are selected as the maxima of the 
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mass distribution within spheres of comoving radius of 1.5h _1 Mpc. The mass of each cluster 
is determined within two relevant radii: a co-moving radius of R com = l-5h _1 Mpc, and a 
physical radius of R p h y = l.Oh -1 Mpc. We use these radii in order to allow a proper compar¬ 
ison with observations, which generally employ R CO m or R p h y as their observable parameter. 
A virial cluster radius, which is commonly used in theoretical analyses (such as the Press- 
Schechter approximation), generally cannot be accurately determined from observations. To 
study the evolution of the cluster mass function, cluster masses are determined at several 
redshifts: z — 0, 0.5, 1 and 2. The cluster MF, n(>M), which represents the number density 
of clusters above a given mass threshold is then determined at each redshift. The evolution 
of the cluster mass function is derived for each model and compared with observations. 

3. Evolution of the Cluster Mass Function 

The evolution of the cluster MF is presented in Fig. 1 for two representative models : 
biased 0=1 SCDM and low-density OCDM, which are “degenerate” at z = 0. A negative 
evolution of the cluster MF is seen in all models — i.e., the abundance of clusters decreases at 
earlier epochs for clusters of a given mass, since massive clusters grow with time (e.g., Press 
and Schechter 1974, Peebles et al. 1989, Peebles 1993, Cen and Ostriker 1994a, Luppino 
and Gioia 1995, Eke et al. 1996). The rate of the evolutionary growth, however, is strongly 
model dependent. For example, the number density of M(< 1.5) > 3 x lO 14 h _1 M 0 clusters 
drops by a factor of ~ 40 from z = 0 to z = 0.5 in the biased SCDM cosmology, while the 
drop is only a factor of 4 (instead of 40) in the low-density models (OCDM, LCDM). The 
difference becomes even larger for more massive clusters (see below). 

The evolution of the cluster density as a function of redshift is presented in Fig. 2 for 
massive clusters with mass M(< R com = 1.5h _1 Mpc) > 5.5 x 1O 14 1 i - 1 M 0 (corresponding 
to richness class ^ 2.5; Bahcall and Cen 1993), and in Fig. 3 for less massive clusters 
M(< R p hy = 1.0h _1 Mpc) > 1.5 x 10 14 h _1 M C | (corresponding to richness class ^ 0, see also 
§4). At z ~ 0, all models except a 8 — 1 SCDM yield a comparable abundance of clusters, 
consistent with observations. The a 8 ~ 1 SCDM model produces an order-of-magnitude 
more clusters than observed. At high redshifts, the abundance of clusters decreases sharply 
for the low a 8 models but the decrease is slow for higher a 8 . The evolution rate is insensitive 
to the value of the Hubble constant, or the exact shape of the power spectrum, and is most 
sensitive to the normalization a 8 (for same mass clusters). The dependence on fl itself is in 
fact only secondary. The strong exponential dependence on <x 8 results from the fact that for a 
given mass cluster, a lower <x 8 implies the clusters are rarer peaks in the density distribution 
, therefore evolving considerably faster than in high <r 8 models (see Fan et al. 1997). An 
observational determination of the cluster evolution rate therefore enables us to directly 
constrain a 8 . 
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4. Comparison with Observations 

Systematic observations of clusters of galaxies at high resdhifts are only now beginning, 
with the use of complete redshift surveys (determining cluster mass from velocity dispersion), 
X-ray observations (temperatures of clusters), and weak gravitational lensing. New complete 
surveys of optical and x-ray clusters at low to high redshifts (z ^ 0.5) will become available 
over the next several years. Here we present results from two independent current optical 
cluster surveys in the redshift range z ~ 0 to ~ 1. While the current samples are still small 
and the uncertainties large, the sensitive cluster evolution already allows us to place strong 
constraints on the cosmology. 

The CNOC optical cluster redshift survey (Carlberg at al. 1996) represents a small 
but complete redshift survey of high mass clusters in the redshift range z = 0.18 — 0.55, 
with an EMSS extension at z — 0.55 — 0.85 (Henry et al. 1992, Luppino and Gioia 1995, 
Carlberg et al. 1997). Redshifts for typically ~ 30 to > 100 galaxies per cluster are used 
to accurately determine the velocity dispersion and mass of each cluster (Carlberg et al. 
1996). The cluster mass threshold and cluster densities in the survey, properly corrected 
for completeness effects, are discussed by Carlberg et al. (1997). The mass threshold used 
is based on a velocity dispersion threshold of a r > 800kms _1 , which corresponds to a mass 
(within R CO m = 1.5h _1 Mpc) of M(< 1.5) > 5.5 x 10 14 h _1 M o (as determined from the 
CNOC data as well as compared with Coma). The cluster densities above this threshold are 
3.5 x 10” 7 Mpc _3 at z — 0.18-0.35, and 9.3 x lO^Mpc" 3 at z = 0.35-0.55 (for Q — 1). The 
high redshift extension at z = 0.55 — 0.85 (Luppino and Gioia 1995, Carlberg et al. 1997) 
does not have complete velocity measurements but contains the richest and most luminous 
X-ray clusters (some with observed velocity dispersions oy ^ 1200km). We conservatively 
assume that these clusters have the same mass threshold as the CNOC clusters (also selected 
from EMSS). The mass threshold is likely to be higher; if so, this will raise the best-fit erg 
value. This uncertainty is included in our estimates. The observed abundance of nearby 
clusters (z ~ 0 — 0.1) is taken from Bahcall and Cen (1993), Mazure et al. (1996, the ESO 
survey) and Henry and Arnaud (1992; based on X-ray cluster selection), all converted to the 
common mass threshold of 5.5 x 10 14 1 i _1 Mq (using the observed MF, Bahcall and Cen 1993). 
This common mass threshold is slightly lower than a Coma-type cluster (~ 6.5 x 10 14 h _1 M Q j 
Hughes 1989, Bahcall & Cen 1993). The sensitivity of the results to the exact mass threshold 
is tested by varying the assumed threshold from 5 to 6.5 xl0 14 h _1 M Q ; the results (§5) include 
these uncertainties. 

The results are presented in Fig. 2, together with the model expectations for this mass 
threshold clusters. (The model a r and mass thresholds are corrected for the resolution effect 
of the simulation by comparing with high resolution simulations; the effect is small: £5%). 
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Only a mild negative evolution is observed. This mild evolution is in excellent agreement 
with the low-density high-normalization models (OCDM, LCDM); it is inconsistent by a 
factor of 10 - 100 with the very strong evolution expected in the biased (<r 8 ~ 0.5) SCDM 
and MDM models. The expected cluster density decreases by a factor of ~ 10 from z = 0 to 
z = 0.5 in the low-density models, while the decrease becomes enormous ( rsj 10 3 ) for a 8 — 0.5 
SCDM, and ~ 10 2 for MDM. The data show a decrease by a factor of ~ 10 to z ~ 0.5. This 
comparison differentiates the D = 1, cr 8 ~ 0.5 models from the D ~ 0.3 — 0.4, <r 8 ~ 0.8 models, 
which are indistinguishable at z ~ 0. Only the low-density, higher normalization models are 
acceptable at high redshifts (see §5). The unbiased a 8 ~ 1 SCDM model, which also yields 
mild evolution (due to its high <r 8 ), is inconsistent with the observed cluster abundance at 
any redshift. 

The second cluster sample we investigate is the Palomar Distant Cluster Survey (PDCS; 
Postman et al. 1996). The PDCS is a complete automated survey of distant clusters to 
z ~ 1 from deep imaging CCD data over 5 deg 2 . Clusters were selected from the imaging 
data using a matched-filter algorithm, which yields best-fit estimates of the cluster richness 
(oc luminosity) and redshift. While the clusters do not have measured redshifts and velocity 
dispersions (i.e., masses), the estimated luminosities are determined in a consistent manner 
from z ~ 0.2 to ~ 1, enabling us to investigate the evolutionary trend of the cluster densities. 
Measurements of the cluster redshifts and velocity dispersions will eventually provide more 
accurate results. We select all clusters with luminosities L d > 50 L* (selected in the /-band, 
with richness threshold A c ; > 50 where L d = A d L* within a physical radius of lh 1 Mpc, 
the radius used by the PDCS selection). This corresponds to a conservative mass threshold 
of M(/? p hy = lh _1 Mpc) ^ 1.5 x 10 14 h _1 Mpc for an average cluster M/L ~ 300h (Bahcall et 
al. 1995, Carlberg et al. 1996). (The evolution results are not sensitive to the exact mass 
threshold at these low mass values; see below). 

Figure 3 presents the evolution of the PDCS cluster density to z ~ 1. The most distant 
point, at z — 0.9, includes an incompleteness correction using the PDCS calibrated selection 
function at that redshift, and is thus less accurate. (The error bar includes both statistical 
uncertainties and a conservatively estimated uncertainty due to the selection correction). 
Figure 3 compares the observed evolution with the model expectations for the same physical 
radius and mass threshold clusters. The data, again, show only a minimal evolution of 
the cluster density, in excellent agreement with the low-density, high-a 8 models. The data 
are inconsistent with the biased SCDM and MDM models, which predict ~ 10 times lower 
cluster density than observed at z ~ 1. (If the actual mass threshold of the clusters is larger 
than estimated above, then the model evolution will be even stronger and the evolutionary 
difference among models somewhat larger; in that respect the assumed mass threshold is a 
conservative choice). 
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The results are similar for the PDCS and the CNOC samples. While the CNOC clusters 
represent considerably higher mass clusters, which are most sensitive to the cosmology, the 
PDCS clusters reach to higher redshifts of z ~ 1. The fact that both independent samples, 
with different mass threshold clusters and different selection algorthiums, yield similar results 
provides further support to these conclusions. 


5. Constraining D 

A comparsion of the observed cluster evolution with the models shows that the data are 
consistent with the low-density models (OCDM and LCDM), and are inconsistent with the 
0=1 models (SCDM, biased SCDM and MDM). The relatively mild evolution observed 
in both the CNOC and PDCS samples is consistent with OCDM at a significance level of 
~ 60% (based on a y 2 test), and with LCDM at ~ 30%. The 0 = 1 SCDM and MDM 
models are rejected at > 99.9%. 

We use the data to directly determine the best-fit values of O and a 8 for the CDM 
models. We use the method described by Fan et al. (1997), correlating the evolution rate 
(n(z)/n(0)) with a 8 to determine cr 8 directly, since the primary dependence of the evolution 
rate is on cr 8 (for same mass clusters). The evolution rate is exponentially dependent on a\ 
- increasing strongly as o 8 decreases; it is nearly independent of other parameters, including 
O (see Fan et al). We determine the best-fit relation for the evolution rate versus <j 8 from 
the model simulations and compare the expected relation with the observed evolution rate. 
We find cr 8 = 0.85 ± 0.15. The observed mild evolution rate of rich clusters thus implies a 
nearly unbiased universe; a strongly biased universe (cr 8 < 0.7) is unlikely since it produces 
considerably stronger evolution than observed. The results are consistent with those of 
Carlberg et al. (1997) of a 8 = 0.75 ±0.1. Combined with the D — a 8 relation for present-day 
cluster adundance for CDM models (from Eke et al. 1996), we find D = 0.3 ± 0.1 (for A 
= 0), and fl = 0.34 ±0.13 ( for A = 1 — O). The results are presented in Fig. 4. The 
figure illustrates the powerful diagnostic of cluster evolution in detaining and Us; it places 
the strongest constraints yet on these parameters. The independent constraint placed by 
cluster dynamics, fl ~ (0.2 ± 0.07)<r^ 1 (assuming linear bias; Bahcall et al. 1995, Carlberg 
et al. 1996, 1997) is also shown in the figure; it is consistent with the above results, yielding 
fl = 0.24 ± 0.1 for a 8 = 0.85 ± 0.15. This fl range provides the overlap of the constraints 
placed by the cluster abundance evolution and cluster dynamics observations. These results 
suggest that we live in a low-density, low-bias universe. Recent observations suggesting a 
minimal negative evolution of the X-ray cluster luminosity function (Henry et al. 1992, 
Castander et al, 1994, Collins et al. 1997, Nichol et al. 1997, Mushotzsky et al. 1997) are 



also consistent with the above findings. 
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Table 1. Model Parameters 
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SCDM 

1.0 
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0.5 
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SCDM 
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0.5 

0.53 

MDM 
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0.0 
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0.60 

OCDM 

0.35 
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0.7 

0.80 

LCDM 
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Figure 1. The evolution of the cluster mass function with redshift, for cluster masses within 
a co-moving radius i? C om = 1.5h _1 Mpc. (For z=l in SCDM <j 8 = 0.53, Press-Schechter 
appoximation is used.) 
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Figure 2. Observed vs. model cluster abundance as a function of redshift for clusters with 
mass M(< R com = 1.5h _1 Mpc) > 5.5 x 10 14 h _1 Mg. The observed abundance at z ~ 0 are 
from Bahcall and Cen (1992), Mazure et al. (1996) and Henry and Arnand (1992). The data 
at z ~ 0.27 and 0.45 are from the CNOC survey (Carlberg et al. 1997), and at z ~ 0.7 from 
Luppino and Gioia (1995). The different symbols represent the observed number densities 
for 0=1 (filled circles), 0=0.35, A=0 (open circles), and 0=0.4, A=0.6 (triangles). 
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Figure 3. Observed vs. model cluster abundance as a function of redshift for clusters with 
mass M(< R phy = 1.0h _1 Mpc) > 1.5 x lO 14 h _1 M 0 . The data are from the PDCS survey 
(Postman et al, 1996). The different symbols for the data are the same as in Fig. 3. 
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Figure 4. Observational constraints on 0 and ct 8 . The band erg = 0.50O _0 ' 50+0 ' 12n represents 
the range due to the present day cluster abundance (for the average of open and A models, 
Eke et al. 1996; the inserted dash-dotted and dotted lines are the best fits for open and A 
models, respectively). The darker band of <t 8 = 0.85 ± 0.15, Q = 0.3 ±0.1 is the const raint 
placed in this paper by cluster evolution (§5). The Q = (0.2 ± O.O?)^ 1 band represents 
cluster dynamics constraint. The dashed lines are the COBE four year data (Bunn and 
White 1996). A low-density low-bias universe with 0 = 0.3 ± 0.1 and ct 8 = 0.85 ± 0.15 best 
fits all the data (darkest region), with cluster evolution providing the tightest constraint. 




